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|  INTRODUCTION 


Background 

Recent  development  of  deployable,  lightweight, 
high-strength  composites  has  generated  significant 
interest  in  their  potential  applications  in  aerospace 
structures,  such  as  the  ground  impact  attenuation 
system  in  the  recovery  technology  of  aerospace 
vehicles  and  in-oibit  construction  of  large  space 
structures1 .  The  composite  consists  of  a  low-density 
fast -deployable  structural  foam  as  the  core  and  a  thin 
fabric/resin  skin. 

Extensive  studies  of  several  composite  systems  of 
this  kind  have  been  conducted  by  the  Air  Force  Flight 
Dynamics  Laboratory  and  the  U.S.  Army  Construc¬ 
tion  Engineering  Research  Laboratory  (CERL)3. 
Many  diverse  aspects  of  the  material  system  have  been 
thoroughly  examined,  ranging  from  the  chemical 
formulation  of  polymer  foams  and  fabrication  proces¬ 
ses  of  composites  to  prototype  testing  of  composite 
structures.  Fundamental  concepts  of  fabrication, 
design,  and  testing  of  the  composite  structure  have 
neen  reviewed.3  The  composite  uses  the  most  pro¬ 
minent  properties  of  the  fabric  material  and  the 
space-filling  foam,  so  that  optimum  performance 
may  be  achieved. 

IT.  Hauler.  "Building  Large  Structures  in  Space.”  Astro¬ 
nautics  and  Aeronautics,  Votume  14  (19761,  pp  55-61;  M.  F. 
Card,  "Trends  in  Aerospace  Structures,”  Astronautics  and 
Aeronautics,  Volume  16  (1978),  pp  82-89;  S  R.  Mehaffie, 
“Foam  Impact  Attenuation  System,"  to  appear  in  A! A  A 
Proceedings  ( 1 979). 

!S.  R.  Mehaffie.  State  of  the  Art  of  Impact  Attenuation 
Concepts  for  RPV  Applications,  Al  l  DL-TR-76-S1  (WPAFB, 
1976);  S.R.  Mehaffie,  Investigation  of  a  Deployable  Poly¬ 
urethane  loam  Ground  Impact  Attenuation  System  for 
Aerospace  Vehicles  (HAS  Tests  A  I-4S).  AFFDL-TR-78-145, 
Volume  I,  (WPAFB,  1977);  S.  R.  Mehaffie,  Investigation  of  a 
Deployable  Polyurethane  Foam  Ground  Impact  Attenuation 
System  for  Aerospace  Vehicles  (FI AS  Tests  #49-91),  AFFDL- 
TR-78-145,  Volume  11  (WPAFB,  1978);  A.  Smith,  Investiga¬ 
tion  of  a  Minimum  Time  Foam  Deployment  System,  draft 
Technical  Report  (C'ERL,  1979). 

3S.R.  Mehaffie,  “Foam  Impact  Attenuation  System,"  A 
Collection  of  Technical  Papers,  AIAA  6th  Aerodynamic  Decel- 
erator  and  Balloon  Technology  Conference,  Houston,  TX, 
March  5-7.  1979. 


Recent  development  of  strong  aramid  fibers  has  led 
to  the  selection  of  Keviar*  fabrics  as  otic  of  (he  most 
ptomising  candidates  lor  the  skin  construction  because 
of  its  excellent  combination  of  mechanical,  thermal, 
and  corrosion-resistant  properties.  Advancement  of 
the  fast-foaming  technology  of  high  polymers  for 
structural  applications  has  made  polyurethane  a  very 
attractive  mate'ia’  for  achieving  the  desired  fabrica¬ 
tion  process  and  structural  performance.  In  fact,  it 
has  been  reported  that,  when  proper  mixing  and 
injection  are  introduced,  polyuiethanc  foams  can  be 
deployed  into  a  desired  structural  geometry  within  a 
few  seconds. 

Since  the  foam-filled  tKn-skin  composite  structural 
members  are  const-ucted  differently  from  those  of 
conventional  monolithic  and  fiber-reinforced  compos¬ 
ite  materials,  some  fundamental  difficulties  have  been 
encountered  in  the  study  of  tire  behavior  of  the  fabric/ 
foam  composites.  These  include  excessively  large 
shear  deformation  of  the  section,  local  failure  of  the 
thin  skin  by  wrinkling  (or  local  buckling),  and  highly 
nonlinear  response  and  creep  behavior  of  the  soft 
polyurethane  foam  core.  These  difficulties  provide  a 
great  challenge  in  the  design  and  analysis  of  the  thin- 
skm,  foam-filled  composite  structures  and  need  to  be 
fully  investigated  for  future  applications  of  the  com¬ 
posite. 

Objective 

The  overall  objectives  of  this  study  were: 

1.  To  reduce  the  operating  time  of  u  deployable 
plastic  foam  system  to  5  seconds,  to  develop  a  mixing/ 
dispensing  system  for  the  foam,  to  evaluate  a  tange  of 
geometric  shapes  obtainable  with  the  foam,  and  to 
evaluate  the  stress/strain  characteristics  of  the  foam 
at  various  times  during  the  foam  curing  process. 

2.  To  evaluate  the  peiformance  of  a  high-strength 
fabric  form  filled  with  polyurethane  foam  to  deter¬ 
mine  experimentally  and  analytically  the  roles  of  each 
component  of  the  composite  material  as  a  potential 
structural  material. 

The  objective  of  this  volume  is  to  present  the 
experimental  and  analytical  study  of  deforma.ion  and 
local  failure  of  polyurethane  foam-filled,  thin  Kevlar/ 
urethane  skin  composite  cylinders  subjected  to  trans¬ 
verse  loading  (four-point  bending). 

Approach 

Chapter  2  describes  the  experimental  procedure.  In 
Chapter  3,  methods  of  analysis  for  studying  linear  and 
nonlinear  behavior  of  the  composite  cylinder  are  pre- 

‘Dupont  register  trademark  for  Aramid  fibers. 
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sented  and  the  local  buckling  criterion  during  com¬ 
posite  deformation  is  derived.  Both  pre-local-buckling 
and  post-local-buckling  responses  of  the  beam  arc 
examined.  Results  obtained  from  the  analyses  arc 
compared  with  experimental  data  and  arc  given  in 
Chapter  4.  Their  implications  on  the  performance  of 
the  composites  arc  discussed.  Effects  of  different 
fabric/urcthane  skins  and  foam  materials  are  studied. 
Emphasis  is  placed  on  determining  the  significance 
of  transverse  shear  deformation  and  nonlinear  core 
response  after  local  buckling  of  the  composite  skin. 


2  EXPERIMENTAL  WORK 

Design  of  the  Experiment 

Cylinders  were  selected  for  study  since  they  repre¬ 
sent  the  simplest  structural  shapes  to  make  without  the 
use  of  elaborate  molds  or  form".  They  can  be  tested 
in  flexure  to  evaluate  the  function  of  each  material 
component.  It  was  recognized,  however,  that  cylindri¬ 
cal  beams  are  not  optimum  shapes  in  a  structural  sense. 

Two  high-strength  fabrics  weie  selected  for  the 
study.  Both  were  woven  of  aramid  fiber  and  coated 
with  polyurethane  elastomer  to  reduce  porosity.  The 
lightweight  fabric  was  Kevlar  29  and  the  heavyweight 
fabric  was  Kevlar  49*.  The  aramids  have  a  very  high 
modulus  of  elasticity  and  high  tensile  strength. 

Two  densities  of  polyurethane  foams  were  selected. 
A  2  Ib/cu  ft  (32  kg/m3)  and  a  4  lb/cu  ft  (64  kg/m3) 
foam  were  to  be  used  in  combination  with  each  fabric 
weight. 

Two  diameters  of  cylinders  were  used  in  the  study. 
The  original  plan  was  to  use  1-  and  2-ft  (0.3-  and 
0,6-m)-diameter  beams  to  vary  the  beam  section. 
When  it  was  found  that. heat  generated  in  the  larger 
beams  by  the  foam  caused  damage,  6-in.  (0.1 5-m)- 
diameter  cylinders  were  used. 

Specimen  Preparation 

The  cylinders  were  made  by  hand  mixing  pre¬ 
weighed  batches  of  foam  materials  and  pouring  them 
into  the  sewn  fabric  forms.  The  expansion  of  the 
loam  filled  the  forms.  The  seams  of  the  forms  were 
arranged  so  that  they  were  along  the  sides  of  the  beams 
jnd  would  not  interfere  with  the  loading.  The  top  of 
each  beam  was  identified  so  that  load  application 
would  be  parallel  to  the  direction  of  rise  of  the  foam. 


♦Kevlar  29  end  Kevlat  49  are  Dupont  registered  trade¬ 
marks. 


The  beams  were  made  without  external  formwork. 
The  fabric  skm  defined  the  beam  shape.  In  the  case 
of  the  I -ft  (0,3-r.i)-diameter  cylinders,  the  beams  were 
very  uniformly  circular  in  cross  section.  The  larger 
beams  were  somewhat  elliptical,  being  about  4  in. 
(102  mm)  less  in  height  than  width.  The  small  (6  in. 
(0. 15m))  beams  were  circular  in  cross  section. 

The  I -ft  (0.3-m)*diameter  beams  were  made  first. 
All  the  necessary  materia!  for  the  fosm  for  each  beam 
was  mixed  in  a  single  butch  and  poured  into  the  form 
at  one  end.  The  opening  in  the  form  for  entry  of  the 
foam  mixture  was  zippered  so  that  it  could  be  securely 
closed  during  the  foaming  process.  A  small  air  vent  was 
provided  in  each  end  to  allow  air  inside  the  form  to 
escape.  Two  beams  of  each  rabric  weight  and  each 
foam  weight  were  made. 

The  2-ft  (0.6-m)-diameter  beams  were  made  in  the 
same  manner  as  the  1-ft  (0.3-m)-diamcter  beams, 
except  half  the  required  material  was  mixed  and 
poured  into  each  end  of  the  fabric  forms.  The  lower 
density  foam  beams  were  made  fust  and  the  elliptical 
cross  section  mentioned  previously  was  noted. 

About  30  minutes  after  the  end  of  foam  expansion 
in  the  higher  density  foam  beams,  the  fust  beam 
exploded.  Or.e  seam  was  split  completely  from  end  to 
end  and  both  circular  end  seams  were  split  about  two- 
thirds  of  the  way  around  (Figure  1 ).  Shortly  before  the 
explosion,  internal  temperature  was  noted  to  be  at 
375°F  (190°C)  and  rising.  The  vast  amount  of  heat 
generated  could  not  escape  and  caused  a  combination 
of  softening  of  the  foam  and  thermal  expansion  of 
both  the  foam  and  the  blowing  ageni  within  the  cells. 
The  use  of  the  2-ft  (0.61-m)-diameter  beams  was 
abandoned  because  of  the  poor  shape  achieved  and  the 
excessive  heat  caused  by  the  section  thickness  and 
quantity  of  foam  being  reacted.  The  ruptured  beam 
showed  extensive  charring  in  the  foam  (Figure  2). 

The  6-in.  (0.15-m)-diameter  beams  were  made 
without  incident. 

Tasting 

The  beams  were  tested  in  flexure  by  a  four-point 
loading  method.  The  reaction  points  were  placed  6  ft 
(1.83  m)  apart  and  the  two  load  Joints  were  evenly 
spaced  2  ft  (0.61  m)  apart  and  2  ft  (0.61  in)  from  the 
reaction  points.  Figure  3  shows  the  loading  arrange¬ 
ment  used  with  all  of  the  test  specimens. 

Load  spreading  saddles  were  fabricated  to  use  in  the 
tests  to  prevent  local  crushing  (Figure  4).  The  load 
bearing  test  set  was  designed  to  simulate  an  antici¬ 
pated  service  loading  case  in  which  two  cylinders 
would  be  in  contact  perpendicular  to  each  other.  The 


saddles  were  3  in.  (76.2  mm)  wide  and  conformed 
to  the  shape  of  the  foam/fabric  cylinders  halfway 
around  their  circumference.  Thee  was  some  crushing 
of  the  foam  by  the  saddles,  especially  in  the  lower 
density  foam,  but  the  effect  was  minimal.  The  effect 
of  the  saddles  on  the  bending  moment  of  the  beams 
was  considered  negligible. 

The  equipment  used  in  the  test  used  a  servoram 
with  a  6-in.  (0.1 5-m)  travel.  Strain  gauges  mounted 
cn  the  fabric  of  the  1-ft  (0.305-m)-diameter  beams 
(Figures  5  and  6)  failed  to  provide  significant  strain 
information.  Ram  stroke  was  recorded  to  provide 
deformation  data.  In  all  except  one  case,  the  beams 
bent  to  allow  full  ram  travel.  Figure  7  shows  typical 
bending  of  a  test  beam.  The  lightweight  fabric  with 
the  higher  density  foam  failed  by  fracture  within  the 
"constant  moment”  portion  of  the  beam  between  the 
load  points. 


3  FORMULATION  OF  ANALYSIS 

The  polyurethane  foam-filled  composite  beam 
studied  in  this  investigation  is  shown  in  Figure  8. 
The  composite  structure  has  a  very  intriguing  combina¬ 
tion  of  physical  properties.  It  is  light,  stiff,  strong, 
provides  excellent  impact  resistance,  and  can  be  shaped 
Figure  4.  Foam  crushing  by  loading  saddle.  efficiently.  The  skin  is  a  thin-layer  composite  consist- 
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ing  of  high-strengtli  Kevlar  fabric  and  polyurethane 
resin.  However,  it  has  a  relatively  low  resistance  to 
compressive  stresses  and  can  lose  its  load-bearing 
capacity  due  to  local  buckling.  The  core  is  rather  soft 
and  deforms  readily  under  stresses.  Since  the  skin 
yields  less  readily  than  the  core  foam,  almost  all  tensile 
and  compressive  forces  initially  are  taken  by  the 
fabric/urethane  skin;  very  little  force  is  by  the  core, 
which  acts  predominantly  in  shear.  As  the  external 
load  increases,  local  buckling  occurs  in  the  compressive 
part  of  the  skin.  The  core  starts  to  take  more  and  more 
compressive  and  shear  stresses  and  deforms  nonlin- 
early. 

The  major  difference  compared  to  the  analysis  and 
design  of  conventional  structures  is  the  need  to 
account  for  the  significant  amount  of  shear  deforma¬ 
tion  and  the  highly  nonlinear  response.  Shear  deforma¬ 
tion  in  conventional  structures  is  generally  very  small 
and  is  usually  neglected,  and  properties  of  the  materials 
are  little  affected  by  sustained  load.  In  the  foam- 
Filled,  thin-skin  structures,  not  only  is  the  shear  defor¬ 
mation  appreciable,  but  the  nonlinearity  of  the  core 
and  the  local  buckling  of  the  skin  show  tremendous 
effects  on  the  load-bearing  capacity  of  the  composite. 
The  analyses  developed  in  this  chapter  are  to  take  these 
effects  into  consideration  and  to  provide  a  means  of 
predicting  the  performance  of  the  composite  beam.  It 
will  be  shown  that  in  the  nonlinear  range  of  deforma¬ 
tion,  the  situation  is  so  complex  that  a  numerical 
technique  is  required.  The  wrinkling  load  in  the  com¬ 
pressive  region  is  derived  to  determine  the  local  failure 
of  the  thin  skin.  Deformation  of  the  composite  beam 
after  local  buckling  is  further  analyzed  using  a  non¬ 
linear  finite  element  method. 

Linear  Analysis  of  Cylindrical 
Composite  Beam  Deformation 

In  the  early  state  of  loading,  small  deformation  in 
the  linear  stress/strain  range  of  the  foam  material  is 
generally  observed.  Consider  a  differential  element,  dz, 
of  the  composite  beam  with  a  circular  cross  section  as 
shown  in  Figure  8.  The  classical  Kirchhoffs  plane- 
remain-plane  hypothesis  is  assumed  in  the  composite 
beam  deformation.  The  axial  strain,  e,  at  any  depth  y 
from  the  central  axis  of  the  beam  can  be  expressed  by 


and 


y*: 


Co 


(Eq  2a] 


0*  =  sin'1  [Eq  2b  | 

Qtn 

where  k*  and  dm  are  defined  as 

dfl  1 

*  K~dz  and  dm  =  2^o  +  di)  [Eq  ?1 

in  which  dj  and  dQ  are  the  inner  and  outer  radius  of 
the  skin  in  the  composite,  respectively. 

The  resultant  axial  force  N  and  the  resultant  bend¬ 
ing  moment  M  acting  on  this  section  can  be  calculated 
by 


N  =  f  (i  E£s)(e„-K*y)  r  dr  d0 

J— 

2 

r2n  -9*  r  — 

+  /  Id2  Et5)(e<»-X*y)  rdr 

Jn-e*  J  — 


and 


d0 


di 

■  f  2  E^ )  (e„ -k  *y)v/dll4y2  dy 
J  v* 


+  |d*  E( r)  (eD  -k  *y jx/d’  -4y 2  dy 
2  [Eq  4] 


=  Ig.  Li  E^5)(e0-K*y)y  rdrd0 
2 

d„ 

f J«+»*  f  f 

+  J„n .  E«5)(eo-K*y)  y  r  dr 


M 


d0 


e  =  e°~Ky  +  dfy 

where  eQ  is  the  axial  strain  caused  by  stretching,  k 
is  the  curvature  of  the  deformed  central  plane,  y  =  0, 
and  0  is  the  angle  of  shear  deformation  of  the  plane. 
The  neutral  axis  and  the  angle  0*  associated  with  the 
axis  as  shown  in  Figure  8  are,  therefore,  located  at 


di 


■  f  2  E£f)(e0-K*y)  y  Vdf-4y2"  dy 


+  /  dj  E|f)(eo-K*y)  y  Vd?-4y2  dy 

"T  [Eq  5] 


) 


t ! 


where  £  is  the  Young’s  modulus,  the  superscripts  (s) 
and  (0  represent  the  modulus  associated  with  the  thin- 
skin  composite  layer  and  tire  foam  core,  respectively, 
and  ;he  subscripts  c  and  t  refer  to  the  modulus  in 
compression  and  tension,  respectively. 

Carrying  out  the  integration  in  Eqs  4  and  5 

N  =  C>  i  c.„  +C|2  x*  [Eq  6] 

MsCn  efl  tCjj  k*  {Eq  7] 

where 

C„  =  J  E^(d„-df )  + 1  (E[f>  +  E<f))  d? 

[Eq  8a] 

c3  2  =  fA  E<s)  (dS  -df  )  +  *“?  (Etf)  +  E(cf) ) 


^(E^-E^df^sin 
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vMF) 

{Eq  8b] 


and 


C,2  =  yjj(E(f)-E£f))(d?-4y*V  [Eq  8c] 

It  should  be  noted  that  bending  and  stretching  are 
coupled  in  the  composite  beam  deformation  if  the 
tensile  and  compressive  moduli  of  the  polyurethane 
foamed  core  are  different.  For  simplicity  without  loss 
of  generality,  the  foam  is  assumed  to  be  isotropic  in 
tension  and  compression;  that  is,  if  Et(f)  =  Et(f) 
=  E*f\  then  the  coupling  effect  disappears. 

C, ,  =5  [E<f>  df +  E(S)  (dj-df)]  [Eq  8a) 

C23  - h  [E(f)  df  +  E<s>  (d*-df )  ]  [Eq  8b] 


and 


C,2  =0 


[Eq  8c] 


Using  the  fundamental  beam  equation,  the  resultant 
shear  force  V  acting  on  a  cross  section  can  be  obtained 
by 


V  = 


dM 

dz 


(Eq  9] 


The  shear  deformation  of  the  composite  beam  may  be 
evaluated  by  considering  that  all  the  shear  force  is 
taken  by  ihe  core,  since  the  cross-sectional  area  of  the 
thin  composite  skin  is  much  smaller  than  that  of  the 
core  (t/d;~0[10-’]). 


0  =  -; 


dM 


[Eq  10] 


G(f)A  dz 

where  and  A  are  the  shear  modulus  and  cross- 
sectional  area  of  the  foam  core,  respectively,  and  k 
is  a  conversion  numerical  factor  having  the  value  of  4/3 
for  a  circular  cross  section  as  suggested  by  Mindlin.4 
Since  tne  linear  deformation  of  the  composite  beam 
w(z)  can  be  related  to  bending  moment  M  by 


M  ~  k*Z 


/d^w_dg\. 
V  dz1  dz  / 


[Eq  11] 


where  Z  is  the  bending  rigidity  of  the  composite  beam 
defined  as 

Z»£[d?  E(f)  +  (d4 -df)E(f)]  [Eq  12] 

for  the  case  of  Eclf*  =  E/f)  =  E(f\  the  deflection  of 
the  beam  can  be  directly  obtained  by  integrating  Eq  1 1 


w(2).//| 


dzdz  - 


16  M 


«  V  HI  _  .  , 

3ffd?G(f)  [Eq  13J 


For  a  four-point  bending  problem  of  a  circular  cross- 
sectional  composite  beam  as  shown  in  Figure  9,  the 
deflection  w(z)  can  be  calculated  exactly  as 


w(z)  = 


PL3 /  z2  23  \  16  PL  _ 

0<z<2 
[Eq  14a] 


PL3 

Z 


1  lz  15  \  16P(2L"^ 


(  z3  +^I—  _i-  11') 

V6L5  4L1  8  L  ~  16/ 


L  3 

-<z<-L 
2  2 


'  3tt  df  G(t> 
[Eq  14b] 


Local  Buckling  (Wrinkling)  of  Composite  Cylinders 

When  the  thin-skin  composite  beam  is  subjected  to 
bending  and  transverse  loading,  nonuniformly  dis¬ 
tributed  axial  compressive  and  tensile  stresses  occur 


4R.  D.  Mindlin,  '‘Influence  of  Rotatory  Inertia  and  Shear 
on  Flexual  Motions  of  Isotropic  Elastic  Plates,”  Journal  of 
Applied  Mechanics,  Volume  18  (1951),  pp  31-38. 
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in  addition  to  transverse  shear  deformation.  The  axial 
compression  and  the  transverse  shear  deformation  in 
the  composite  may  introduce  local  buckling  (or  wrinkl¬ 
ing)  of  the  composite  skin  in  which  the  transverse 
normal  stiffness  may  play  an  essential  part.  This  local 
buckling  phenomenon  is  a  form  of  instability  assoc¬ 
iated  with  shoit  waves  of  skin  deformation,  which 
under  a  compressive  load  may  be  cither  symmetrical 
(out  of  phase)  or  aniisymmetrical  (in  phase)  with 
respect  to  the  center  plane  of  the  core.  The  wrinkling 
is  very  pronounced  and  is  of  special  importance  for 
sandwich  composite  construction  having  a  continuous 
low-density  core  material  such  as  the  soft  polyurethane 
foam.  Local  buckling  of  the  beam  results  in  the  loss  of 
load-bearing  capacity  of  a  part  of  the  thin  composite 
skin  under  compression,  while  the  overall  structure  is 
still  capable  of  taking  external  loads.  The  load-carrying 
cross-sectional  area  of  the  composite  member  progress¬ 
ively  reduces  as  the  wrinkling  in  the  Kevlar  fabric/ 
urethane  skin  develops  when  external  loading  increases 
continuously. 

For  a  very  thin-skin  cylindrical  composite  beam 
(t/dj  <  50),  the  approximate  method  originally 
suggested  in  Williams5  is  applicable  to  the  present 
problem.  By  evaluating  the  total  potential  energy  of 
the  system  (i.e.,  bending  energy  of  the  skin,  strain 
energy  of  the  core  associated  with  the  normal  and 
transverse  shear  stresses,  and  the  potential  energy  of 
the  external  load)  and  using  the  standard  minimiza¬ 
tion  procedure  for  an  elastic  stability  problem,  the 
critical  wrinkling  r, tress  cb  can  be  obtained  by6 

ob  3-0.76  [G(f)E(f)E(s)] 3  [Eq  15] 


Local  buckling  will  occur  and  its  effects  should  be 
considered  in  the  deformation  of  the  composite 
beam  when  maximum  compressive  stress  in  the  beam, 
corresponding  to  an  increasing  external  loading  Pb, 
reaches  this  value. 

Consider  the  present  four-point  bending  problem 
of  the  composite  beam  as  shown  in  Figure  9.  Maximum 


5D.  Williams,  Sandwich  Construction-A  Practical 
Approach  for  the  Use  of  Designers,  RAE  Report  No.  Struc¬ 
tures,  2,  M.  2466(1947). 

6 A.  C.  Eringen,  "Buckling  of  a  Sandwich  Cylinder  under 
Uniform  Axial  Compressive  Loading,”  Journal  of  Applied 
Mechanics,  Volume  18,  No.  2  (1951),  pp  195-202. 


bending  occurs  in  the  middle  section  of  the  span 
between  z  =  -L/2  and  z  =  L/2.  The  maximum  com¬ 
pressive  stress  omax  occurs  at  the  point  (x  =  0,  y 

=  °f  any  section  in  the  middle  span  and  has  the 

value 


°max 


=  E(5>  e0-(K- 


d&do 
dz  2 . 


[Eq  16] 


The  pure  bending  condition  in  the  midspan  of  the 
cylindrical  beam  implies 


and 


[Eq  17] 


[Eq  18] 


Thus,  the  critical  external  force  Pb  at  which  local 
buckling  of  the  composite  skin  occurs  can  be  shown 
to  have  the  following  form 


Pb~iApj[G(f)E(f)E(S>1'  lEql9] 

Post-Local-Buckling  Behavior 
of  Cylindrical  Composite  Seams 

When  the  composite  beam  is  subjected  to  a  four- 
point  bending  load,  local  buckling  generally  starts  in 

the  upper  surface  (x  =  0,  y  =  of  the  filled  cylinder, 

where  the  compressive  stress  is  maximum  and  reaches 
the  local  buckling  stress  ob  as  discussed  in  the  pre¬ 
vious  section.  As  the  external  load  increases,  a  larger 
portion  of  the  composite  skin  undergoes  compressive 
stress  exceeding  ob ;  thus,  the  wrinkled  area  progress¬ 
ively  inceas^s.  While  local  buckling  in  a  composite 
skin  continuously  extends  to  a  larger  dimension,  the 
nonwrinkled  part  of  the  section  can  still  take  further 
increasing  loads.  In  the  analysis  of  post-buckling 
deformation  of  the  composite  beam,  it  is  convenient 
to  divide  the  deformed  cross  section  area  into  three 
distinct  regions  (Figure  10):  (1)  wrinkled  portion  in 
the  skin  (y  >  yb),  where  load-bearing  capacity  of  the 
skin  has  been  lost  entirely;  (2)  the  active  compressive 
region  (yb  >  y  >  y*),  where  both  the  skin  and  the 
core  are  subjected  to  a  compressive  stress  a  (|o|  < 
|ob|)  with  Ec(s)  =  E(s)  and  Ec<f)  =  E(,);and  (3)  the 
region  under  tensile  stresses  (y  <  yb)  where  E,*s*  = 
E^8)  and  E/f)  =  E*f^.  The  lesultant  axial  force  N  and 
the  resultant  bending  moment  M  can  be  evaluated 
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similarly  to  those  in  Eqs  4  and  5,  except  the  first  term 
in  the  right-hand  side  of  Eq  4  is  replaced  by 


fe b  /-do/2 

2/  /  E(,)  (e„  -  x*y)  r  dr 

JB*  4>|/2 


dfl 


=  ^(flb-fl*)(d.^-df)e0 

C(»)  ,  , 

+  -  d?)(cosflb  -cos0*)k* 

[Eq  20] 


where 


and 


6b  =  sin 


in-i  f _ 

l(d,  +  d0 


)/2 


yb  = 


- 


Ob 


°  “  E(s) 


Ik* 


[Eq  21a] 


[Eq  21b] 


and  the  first  term  in  the  right-hand  side  of  Eq  5  is 
replaced  by 


ndo/2 

E<s>fe0-*c*y)y  rdr 
i/a 


dfl 


E<s) 

=  -jy  (do  -  df )  (cosflb  -  cos0*)eo 

p(s)  1 

+  64  (d"  "  d‘  )[Gb  "  **  -  f(sin2flb  -  sin20*)]tc* 

[Eq  22] 

Thus,  one  obtains  the  following  expressions: 

N  =  Cu  e0  +  Cu  k*  [Eq  23a] 

M  =  C2i  e0  +C„  k*  [Eq  23b] 

where 


p(») 

c»i  =  §“(do  -  df) 


*  +  2sin' 1 

dm  J 


[Eq  24a] 


c(»)  i 

c>2  =c2i  =-j|y-(dJ-d?)(dJm-4y6)2 


Cl2  “ 


trE(f> 
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E(l) 

+  54"  (do  ~  df  ) 


i 

n  .  ,  2yb  2yh(d*1-4yg)J 

r  +  sin  - - 7j - 

2  dm  dm 


[Eq  24c] 

Solving  e„  and  x*  from  Eqs  23a  and  23b,  the  post- 
buckling  deflection  w(x)  of  the  cylindrical  composite 
beam  can  be  calculated  by 


w(x)=  [  l\*Uy)  d*,d$ 
j  0  Jo 


-16  PL 


3ffd?G(f) 

-  1 6  P  (|  L  -  z) 
3  ix  d?G<f> 


0<z«^ 


L  3 
2  < z  <- L 


[Eq  25] 


It  should  be  noted  that  nonlinear  deformation  of  the 
soft  core  usually  occurs  after  local  buckling  of  the 
composite  skin.  Thus,  the  nonlinear  constitutive  rela¬ 
tion  of  the  core  material  should  be  included  in  the 
post-buckling  analysis.  A  highly  nonlinear  governing 
differential  equation  is  therefore  obtained.  It  can  be 
solved  only  by  the  numerical  methods  discussed  in  the 
next  section. 

Nonlinear  Analysis  of  Cylindrical 
Composite  Beam  Deformation 

The  rigidity  of  polyurethane  foam  in  the  core  is 
much  smaller  than  that  of  the  load-bearing  Kevlar 
skin  in  the  deployable  composite  structure.  Core 
deformation  in  the  foam-filled  cylinder  under  trans¬ 
verse  loading  is  observed  generally  deviating  significant¬ 
ly  from  linearity,  especially  after  local  wrinkling  of 
the  skin  starts.  To  understand  the  nonlinear  behavior 
better  and  to  obtain  realistic  quantitative  results  for 
design  and  comparison  purposes  (with  experiments), 
it  is  necessary  to  include  the  post-yielding  behavior  of 
the  foam  material  in  the  study  and  to  conduct  a  non¬ 
linear  analysis  of  the  problem.  In  the  nonlinear  range 
of  composite  deformation,  governing  differential 
equations  become  nonlinear  and  considerably  more 
complicated  than  those  of  the  linear  case.  Numerica. 
methods  such  as  the  finite  element  analysis  are  pre  ¬ 
ferable.  Other  complications,  such  as  local  failure  of 
the  skin  by  wrinkling  and  shear  deformation  can  also 
be  incorporated  easily  into  the  formulation. 
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[Eq  24b] 


Consider  material  nonlinearity  of  the  foam  by 
expressing  its  constitutive  equation  in  a  power  law 
form  as  core  deformation  exceeds  the  yield  strain 


rr<f>  = 


f  E(f)  e 
A(e  -  0.02)n 
— A(|e  +  0.02|)n 


kK  0.02 
e  >  0.02 
e  <  -0.02 


[Eq  26] 


where  the  coefficient  A  and  the  hardening  exponent  n 
can  be  determined  from  experiments.  The  Kevlar 
composite  skin  is  assumed  to  deform  linearly  in  both 
tension  and  compression  until  failure.  Thus,  the  in¬ 
cremental  skin  and  core  stresses  6o<s*  and  5o^f\  can 
be  written  as 


8o(s*  =  E(s)  (5e„  -  y  5  k*)  [Eq  27] 


and 


6  a 


(0  = 


E<f) (Se„-yS  «*) 

nA(e-0.02)nl  (5en -y  8  k*)  y,  <y<- 


yt  <  y  <  yt 

di 

2 


-nA(|e  +  0.02|)nl  (8ec -ySs*)  ^  <  y  <  yc 


[Eq  28] 

where  yc  and  yt  are  planes  in  a  section  at  which  local 
compressive  or  tensile  strains  reach  its  yield  strain 
(eys  =  ±0.02).  The  corresponding  incremental  resultant 
force  8N,  bending  moment  6M.  and  shear  force  6v  can 
be  obtained  by 


6N=  j 

i  80  dA 

[Eq  29a] 

J  A 

6M=  1 

f  -  y  da  dA 

[Eq  29b] 

J  A 

6V  =  |g<0  A  6/3 

[Eq  29c] 

and 


Incremental  Element  Stiffness  Matrix  Formulation 

In  the  formulation  of  the  nonlinear  finite  element 
analysis  of  composite  beam  deformation-consider  a 
cylindrical  beam  element  (Figure  11)  subjected  to 
incremental  loading-the  stress/strain  relation  can  be 
written  as 


8N 

c„  Cu  0 

6e„ 

8M 

= 

C21  c22  0 

6k* 

5V 

0  0  C33 

6/3 

[Eq  30a] 


or 


where 


6a  =  C  6e 


[Eq  30b] 


Ctl  =gE(s>  (d,2,-d2)^jr  +  2  sin'1  —) 

/•y  c  _ 

+  /  E<f)  Vdf-4yJ  dy 
•'yt 

/-yt  , _ 

+  I  d|  nA(e-0.02)n  l  Vdply2  dy 


/  :  nA(|e  +  0.02|)n  1  Vd?-4yJ  dy 

>c 


[Eq  31a] 


c.2 =  C2.  “  E(s)  (d2-d?)  v^yT 


r  yt  , _ 

I  d.nA(e-0.02)n  ly  Vdf^4y2  dy 


2 

«li 


+  /  2  nA(|e  +  0.02|)n  ly  x/J?-4y2  dy 

yc 

[Eq  31b] 


C22  =^E<‘>  (do— df )  (^  +  sin' 1 


+  f  E(f)  y2  >/df-4y2  dy 
yyt 


and 


r  yt  _ _ 

+  J  dj  nA(e-0.02)n'‘  \/d2-4y2  y2  dy 

~  T 

fT 

-  /  2  nA(|e  +  0.021)"-'  Vd?-4y2  y2  dy 

yt 

[Eq  31c] 

C33  =|c(f)  A  =  -^7rdf  G(f)  [Eq  3 Id] 
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Figure  10.  A  cross  section  of  locally  buckled  composite  cylinder  element. 


The  strain-displacement  equation,  by  definition,  has 
the  following  form 


potation  functions  (linear  for  u  and  0,  and  cubic  for 
Figure  11.  An  element  of  deployable  composite  w),  the  displacement  u  in  an  element  may  be  expressed 

cylinder  for  nonlinear  FEM  analysis.  by  element  nodal  displacements  q 


l-tl  0  0 

o  ti-.u1  +  :tJ>  /(i-t1 
o  o  o 


Of  o  oo 
o  o  f'u-:t)  o 
d-ft  oo  Of 


[Eq  34] 
or 


u  =  Rq  1  Eq  35] 


where  £  =  z/K;  V  is  the  length  of  the  element;  and  ui( 
Wj,  (,  0,,  (i  =  1,2)  are  defined  in  Figure  11.  Thus,  one 
obtains 


incremental  and  iterative  procedures  for  solving  the 
material  nonlinear  problem.  The  overall  procedure 
is  shown  graphically  in  Figures  12  and  13.  The  load 
is  divided  into  many  increments  or  panial  loads  and  is 
applied  one  increment  at  a  time  as  shown  in  Figure 
12.  Within  each  load  increment,  a  tangent-stiffness 
iterative  procedure®  (Figure  13)  is  useu  to  determine 
the  true  solution. 

The  final  solution  after  n  steps  of  increments  is 
equal  to  the  summation  of  all  incremental  solutions. 
This  mixed-mode  solution  scheme  can  be  described 
mathematically  as  the  following.  The  load  Q|  after 
the  application  of  the  ith  increment,  is  given  by~ 


e  =  Bq  [Eq  36] 

where 

B  »  DR  [Eq  37] 

Following  the  routine  procedure  of  the  variational 
principle  of  minimum  potential  energy,  the  element 
stiffness  matrix  k  during  an  incremental  loading  step 
can  be  shown7  to  have  the  form 


and  the  loading  vector  Q  is  defined  as 

QT  =  [F, , P, ,  M, , -M| ,  Fj,  P-. ,  Vj, -M»] 

lEq  39] 

where  superscript  T  represents  the  transpose  of  the 
associated  matrix,  since  the  wmk  done  due  to  external 
loads  is 

W-Fiu,  +  P|  W|  +  M|(0  -  0, )  +  T;u,  +  PjWj 

+  M2(O2-02)  [Eq  40] 

Eqs  38  and  39  _re  used  directly  in  a  mixed-mode 
solution  procedure  for  the  nonlinear  analysis  of  the 
composite  cylinder. 


Qi  =  Qo  +  ji  AQj  [Eq  41] 

An  increment  to  the  displacements  qi  after  nth 
itei  ction  can  be  computed  by 

K}">Aq|n>=(Ql+l-Q,)-t;n-,) 

[Eq  42] 

where  is  the  load  equilibrated  after  the  pre¬ 

vious  cycle,  and  K((n\  the  stiffness  of  the  nth  cycle,  is 
usually  chosen  by  the  tangent  stiffness  matrix  at  the 
end  of  the  previous,  i.e.,  the  (n-l)th  iterative  step. 
Thus,  Eq  42  becomes 

K<n->>  Aqf">=(Qitl  -Qi)-iJn  ° 

[Eq  42] 

Total  displacements  after  nth  iteration  during  a  load 
increment  (Qi+ ;  -  Q,)  are  computed  from 

qjn)=q{0>+  l  Aq[m>  ]Eq  44] 

~  m  “  1 


The  current  equilibrated  load  I  (<n)  is  evaluated  as 
the  load  necessary  to  maintain  the  displacement  qj(n) 


Step-Iteration  Solution  Scheme  for  Nonlinear  Analysis 

The  present  solution  scheme  uses  a  combination  of 


7C.  S.  Desai  and  J.  F.  Abel,  Introduction  to  the  Finite 
Element  Method  (Van  Nostrand  Reinhold  Company,  1977). 


This  iterative  procedure  is  repeated  until  the  incre¬ 
ments  of  displacements  or  unbalanced  loads  become 
zero,  i.e.,  Aq  t<n)  or  (Qj+i-Qi  -Tjn)  sufficiently 


®C.  S.  Desai. 
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close  to  null  according  to  some  preset  criterion.  The 
total  displacements  after  ith  increments  of  loads  are 
calculated  by 

q{  =  4„  +  £  qj  lEq46] 

~  ~  It  - 

Then  additional  incremental  load  (Qi+:  -  Qi+i)  «•' 
added  and  the  similar  procedure  is  repeat*  1  until  the 
applied  load  level  is  reached. 


MECHANICAL  PROPERTIES 
4  OF  FABRIC/URETHANE  COMPOSITE 
SKIN  AND  POLYURETHANE  FOAMS 


High  modulus,  high-strength  Kevlar  aramid  fibers 
offer  an  extremely  desirable  performance  which  is 
well  suited  to  the  construction  of  a  wide  range  of 
engineering  components  and  structures  using  advanced 
liber-reinforced  composites.  The  use  of  Kevlar  fibers 
in  the  form  of  fabrics  and  fabric  composites  has  been 
relatively  new.  Fabrics  of  the  Kevlar  aramid  have  been 
found  to  possess  similarly  attractive  properties  such  as 
high  strength;  low  weight;  good  dimensional  stability; 
excellent  tearing,  cutting,  and  fatigue  resistance; 
thermal  stability;  and  flame  resistance.  Thus,  the 
fabrics  are  uniquely  qualified  for  a  variety  of  special 
load-bearing  elements  in  engineering  applications  of 
the  material  such  as  ballistic  protection  materials, 
inflatable  structures,  and  the  impact  attenuation 
system  in  aerospace  vehicle-recovery  technology.  The 
strong  skin  of  the  fabric/foam  composite  structure  is 
made  of  a  thin  layer  composite  consisting  of  Kevlar 
fabric  and  polyurethane  resin.  Two  kinds  of  commonly 
employed  aramid  fibers,  Kevlar  29  and  Kevlar  49,  are 
presently  used  in  the  skin  construction.  The  fabrics 
have  a  plain  weave  structure  with  a  weight  of  7.0 
oz/sq  yd  (0.24  kg/mJ)  and  10.0  oz/sq  yd  (0.34  kg/mJ), 
respectively.  Fiber  volume  fractions  along  the  warp  and 
filled  directions  in  the  composite  skin  are  both  about 
15  percent.  Young’s  moduli  of  the  composite  skin 
have  been  determined,  both  expeumentally  and 
theoretically,9  giving  values  of  about  1,35  X  106  psi 
(9.31  X  107  Pa)  and  2.85  X  106  psi  (1.96  X  1010  Pa) 


9M.  W.  Waidle,  Textile  Research  Laboratory,  E.  I.  Dupont 
lie  Nemours  &  Company,  Wilmington,  DE,  private  communi¬ 
cation  (1979). 


for  the  Kevlar  29  and  Kevlar  49  fabric  composite  skins, 
respectively.  These  material  properties  will  be  used 
directly  in  the  present  analysis. 

i.i  contrast  tc  the  ductile  nature  of  certain 
polymers,  the  poly  methane  foam  core  tends  to  exhibit 
brittle  behavior  and  has  very  low  strength  in  tensile. 
However,  in  compression,  the  foam  has  constitutive 
behavior  similar  to  ductile  polymers  with  an  apparent 
yield  stress,  large  elongation,  and  relatively  high 
ultimate  tensile  strength.  The  apparent  yield  stress  is 
considered  generally  to  be  related  to  the  collapse  of 
the  cell  structure  in  the  foam.  Although  several 
theories  for  evaluating  the  elastic  modulus  of  the  foam 
have  been  proposed,  the  Kcrncr  equation  and  the 
modified  Haipin-Tsai  equations10  seem  to  give  the 
most  favorable  results  when  compared  to  the  experi¬ 
ments.  For  the  low-density  foams  used  in  this  study, 
elastic  moduli  of  the  foams  have  Ken  found  to  be 
related  to  the  density  of  the  foam  by  a  linear  relation¬ 
ship. 

Figure  14  shows  the  experimentally  determined, 
static  stress/strain  behavior  of  two  kinds  of  polyure¬ 
thane  foams  with  different  densities  (2  and  4  lb/cu 
ft  [32  and  64  kg/m3  J ).  The  nonlinearity  in  the  figure 
is  quite  obvious  as  strains  exceed  2  percent.  For  the 
convenience  of  the  analytical  development,  the  experi¬ 
mental  data  were  fitted  into  continuous  curves  by  a 
best  least-square,  curve-fitting  technique  to  give  con¬ 
stitutive  equations.  In  the  elastic  range  of  foam  defor¬ 
mation,  Young’s  moduli  of  the  foams  have  been  deter¬ 
mined  to  he  850  psi  and  2400  psi  (5.86  X  I06  Pa  and 
1.65  X  107  Pa),  respectively,  for  the  low  and  high 
density  foams.  In  the  inelastic  range,  the  constitutive 
relations  of  the  foams  have  the  following  forms: 

a  =  [48.  +  260.  (|e|-0.02)°  7  J  sgn(e) 

[Eq  47a] 

for  the  4  l^/cu  ft  (64  kg/m3)  polyurethane  foam,  and 
[17.  +  47.2  (|e|-0.02)°  42  J  sgn(e) 

[Eq  47b] 


10G.  J.  Fallick,  H.  J.  Bixler,  R.  A.  Marsella,  V.  F.  Garner, 
and  E.  M.  Fettes,  Modem  Plastics,  Volume  45,  No.  5  (1968), 
pp  143-154;  L.  E.  Nielsen,  Mechanical  Properties  of  Polymers 
and  Composites,  Volume  2  (Marcel-Dekker  Inc.,  New  York, 
1974). 


23 
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Figure  14.  Stress/strain  behavior  of  polyurethane  foams  with  two  different  densities. 


for  the  2  lb/  -i  It  (.12  kghn')  foam.  It  should  be  noted 
that  foams  are  generally  anisotropic  because  gravity 
and  restraint  effects  during  the  foaming  process  tend  to 
elongate  the  foam  cells  in  the  transverse  direction. 
Thus,  it  is  more  difficult  to  collapse  the  cell  structure 
in  the  thickness  direction  than  in  the  longitudinal 
direction.  The  curve  extends  to  exhibit  higher  elonga¬ 
tion,  and  then  the  stress  rapidly  increases  after  the 
cell  wall  has  collapsed,  which  is  not  shown  in  the 
figure. 


5  RESULTS  AND  DISCUSSION 


The  foam-filled  fabric/urethane  skin  composite 
cylinder  of  length  3L(L  =  2  ft  (.607  m| )  and  diameter 
dm2r(dm  =  6  and  12  in.  (1 52  and  305  mm])  subjected 
to  four-point  bending,  as  shown  in  Figure  9,  is  con¬ 
sidered.  Two  kinds  of  fabric  material,  i.e.,  Kevlar  29 
and  Kevlar  4l>,  are  used  in  the  skin.  The  Kevlar  fabric/ 
urethane  skin  has  a  thickness  of  about  0.012  to  0.015 
in.  (0.3048  to  0.3810  mm).  The  mechanical  properties 
are  discussed  in  the  previous  section.  Two  kinds  of 
low  density  polyurethane  foam  (2  and  4  Ib/cu  tt 
(32  and  64  kg/m3])  are  used  as  the  core  material  in 
the  composite  construction.  The  foams  have  very  low 
elastic  moduli  and  deform  rapidly  into  nonlinear  range 
(Figure  14).  It  is  found  that  the  Kevlar  fabric/urethane 
skin  bonds  very  well  with  the  core  material.  Thus,  for 
each  composite  beam  geometry,  four  different  cases 
are  to  be  studied.  In  a  four-point  bending  test,  the  load 
P  increases  continuously,  and  deflection  of  the  com¬ 
posite  cylinder  is  recorded.  Initial  elastic  deformation 
of  the  beam  and  subsequent  local  buckling  (wrinkling) 
of  the  skin  are  examined  by  the  methods  of  analysis 
described  in  Chapter  3.  Progressive  extension  of  the 
wrinkling  portion  in  a  beam  section  until  total  failure 
of  the  composite  by  tensile  fracture  of  the  skin  is 
determined  quantitatively  by  the  nonlinear  step- 
iterative  solution  scheme  mentioned  earlier.  Analytical 
solutions  for  each  case  are  compared  with  experi¬ 
mental  results  to  ensure  the  validity  of  the  analyses  and 
to  assess  the  accuracy  of  the  results. 

Accuracy  and  Convergence  of  Solutions 

Accuracy  and  convergence  of  present  solutions 
obtained  from  the  aforementioned  nonlinear  analyses 
have  been  examined  by  solving  problems  whose  solu¬ 
tions  are  ava'lable  in  the  literature  and  by  comparing 
present  analytical  results  with  experimental  data.  The 


nonlinear  deformation  of  a  cantilever  beam  (Figure  15) 
for  which  plastic  behavior  of  von  Mises  type  is 
assumed,  is  chosen  here.  This  problem  is  attacked  by 
the  present  step-iterative  solution  scheme.  Figure  16 
shows  the  currently  calculated  load-deflection  curve 
and  the  result  obtained  from  Zienkiewicz.1 1  As  is 
expected,  as  the  collapse  load  is  approached,  progress¬ 
ively  large  numbers  of  iterations  are  required.  The  two 
results  appear  in  excellent  agreement  with  each  other. 
The  convergence  of  solutions  for  nonlinear  problems 
depends  on  the  selection  of  load  increments  and  the 
degree  of  discretization  (i.e.,  number  of  elements). 
Figure  1 7  illustrates  graphically  the  rate  of  convergence 
of  present  solutions  by  plotting  the  number  of  ele¬ 
ments  vs.  the  deflection  and  the  curvature  at  the  end 
of  the  center  line  of  the  deformed  beam.  It  :s  observed 
that  an  error  of  about  O.i  percent  in  the  solution  can 
be  achieved  by  using  1 8  elements  or  more  in  the  pre¬ 
sent  analysis.  This  example  gives  important  informa¬ 
tion  concerning  the  discretization  of  the  continuum  in 
the  present  study.  All  subsequent  analyses  for  various 
kinds  of  foam-filled  composite  cylinders  use  numbers 
ot  elements  which  will  guarantee  an  accurate  solution 
within  0.1  percent. 

Prediction  of  Deformation  and  Local  Buckling 
of  Various  Fabric/Urethane  Foamed 
Composite  Cylinders 

The  currently  established  methods  are  now  used  to 
examine  actual  deformation  of  the  composite  cylinders 
under  four-point  bending.  Sohrions  have  been 
obtained  for  the  cases  of  various  combinations  of 
materials  and  geometrical  variables.  Figures  18,  24, 
29,  and  35*  give  both  predicted  and  experimentally 
observed  load-deflection  curves  of  the  1 2-in.  (0.3048- 
m)-diameter  composite  beams  containing  different 
combinations  of  fabric/urethane  skins  and  polyure¬ 
thane  foams.  Figures  37,  39,  43,  and  47  provide 
similar  information  for  the  6-in.  (T52-mm)-diameter 
composite  cylinders  with  various  material  properties 
and  geometric  variables.  Open  circles  in  the  figures 
represent  experimental  data  directly  determined  from 
the  tests.  Solid  lines  indicate  the  response  of  the  com¬ 
posite  beams  predicted  by  the  present  nonlinear 
analysis  with  transverse  shear  deformation  and  local 


11 0.  C.  Zienkiewicz,  The  Finite  Element  Method  in 
Engineering  Science  (McO-aw-Hill,  1971). 

•The  appearance  of  nonsequential  figure  numbers  in  the 
text  is  intentional.  The  grouping  of  basic  data  figures  followed 
immediately  by  supporting  analytical  figures  accounts  for  the 
arrangement. 
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Figure  16.  Deflection  at  the  end  of  the  beam  vs  increasing  load. 


Figure  17.  Accuracy  and  convergence  of  the  nonlinear  FEM  analysis. 


buckling  of  their  skin  included  in  the  formulation. 
These  figures  have  shown  that  the  analytical  predic¬ 
tions  provide  accurate  results  for  almost  all  cases  as 
compared  with  the  experimental  data  when  the  non¬ 
linear  behavior  of  the  core,  the  shear  deformation  of 
the  soft  foams,  and  the  local  buckling  of  the  Kevlar 
fabric/urethane  skin  are  properly  considered  in  the 
analysis.  The  local  buckling  of  the  thin  composite 
skin  starts  at  a  rather  low  load  level  and  has  also  been 
determined  by  the  analysis.  It  is  observed  that  the 
wrinkling  load  Pb  is  related  to  the  elastic  moduli  of 
both  the  core  and  the  skin  materials  and  that  the  core 
material  has  a  greater  effect  on  the  local  buckling 
than  that  of  the  composite  skin.  For  example,  the 
wrinkling  load  for  the  composite  beam  of  Kevlar  49 
fabric/4  lb/cu  ft  (64  kg/m3)  urethane  foam  is  nearly 
twice  as  much  as  that  of  Kevlar  49  fabric/2  lb/cu  (32 
kg/m3)  urethane  fcam  composite. 

Effects  of  Shear  Deformation 

One  of  the  most  important  problems  investigated  in 
the  present  study  is  the  significance  of  shear  deforma¬ 
tion  and  its  effects  on  the  overall  deflection  of  the 
composite  cylinder.  In  a  four-point  bending  test,  total 
deflection  of  the  beam  may  be  conceived  as  composed 
of  three  components,  caused  by  axial  stretching, 
bending,  and  shear  strains.  Shear  deformation,  usually 
considered  negligibly  small  in  conventional  structures, 
has  been  observed  to  be  very  appreciable  in  the  foam- 
filled  composite  subjected  to  transverse  loading.13  This 
phenomenon  becomes  increasingly  pronounced  after 
the  fabric/urethane  skin  starts  wrinkling.  It  is 
extremely  significant  in  the  later  stage  of  the  beam 
deformation,  when  a  large  portion  of  the  core  section 
is  strained  into  the  nonlinear  range.  Present  results  have 
shown  that  shear  deformation  in  the  foam-filled 
thin  skin  composite  plays  an  extremely  important  role 
in  governing  the  overall  response  of  the  composite 
structure.  The  shear  effects  have  to  be  considered  in 
the  formulation  and  analysis  of  composite  cylinders. 
Failure  to  include  the  shear  deformation  in  the  analysis 
will  result  in  a  significant  error  in  the  final  results. 
Figures  IS,  24,  and  29,  for  example,  give  quantitative 
information  of  the  effects  of  shear  deformation.  It 
has  been  found  that  shear  deformation  contributes  to 
about  half  of  the  total  deflection  of  the  cylindrical 
composite  beams  during  the  entire  loading  history. 
If  not  properly  taken  into  account  in  the  study,  the 


,2F.  H.  Wittmann,  Stevin  Report  10-76-5,  Stevin  Labora¬ 
tory  for  Mechanics,  Structures  and  Materials,  Delft  University 
of  Technology,  The  Netherlands,  1976. 


predicted  load-deflection  curves  will  deviate  from  the 
actual  experimental  results  by  a  significantly  large 
amount.  In  fact,  it  has  been  observed  that  if  the  shear 
deformation  is  not  included,  it  will  not  be  possible 
to  predict  the  actual  response  of  the  composite  for 
any  of  the  cases,  regardless  of  whatever  other  com¬ 
plications  are  included  in  the  formulation. 

Core  Stresses  and  Nonlinear  Effects 

When  stresses  in  the  composite  cylinder  are  built 
up  during  loading,  especially  after  local  buckling  of  the 
skin,  Lhe  core  material  starts  to  pick  up  more  loads 
under  increasing  external  loads.  Nonlinear  deformation 
in  the  foamed  core  section  gradually  becomes  very 
significant.  The  development  of  core  stresses  can  be 
calculated  precisely  by  the  present  nonlinear  solution 
scheme.  Figures  19  through  23,  26  through  28,  30 
through  34,  40,  41,  42,  44,  45,  and  46  present  quanti¬ 
tative  information  of  the  increasing  core  stresses  from 
relatively  small  and  linear  distributions  in  the  initial 
stage  to  highly  nonlinear  ones  in  the  later  stage  of 
deformation.  The  nonlinear  behavior  of  the  core 
affects  the  total  deformation  of  the  cylindrical  com¬ 
posite  beam,  as  can  be  illustrated  in  the  predicted 
load-deflection  curves  for  various  cases  (Figures  18, 
24,  35,  37,  38,  43,  and  47.  It  is  observed  that  devia¬ 
tions  of  deflection  from  linearity  become  increasingly 
significant  so  that  it  would  be  impossible  to  obtain  an 
accurate  prediction  of  deformation  if  the  nonlinear 
effects  of  the  core  are  not  included  in  the  analysis. 

Progressive  Extension  of  Wrinkling  and 
Compressive  Zones  in  Composite  Sections 

It  has  been  discussed  in  the  previous  section  that 
when  local  buckling  of  the  thin  Kevlar  fabric/urethane 
skin  starts,  the  compressive  load-bearing  area  in  the 
skin  begins  to  decrease  and  the  compressive  zone  in 
the  core  begins  to  increase  so  that  the  equilibrium 
condition  still  can  hold.  As  loading  continuously 
increases,  progressive  extension  of  the  wrinkling  and 
compression  zones  continues.  The  reduction  of  the 
load-bearing  domain  in  the  strong  skin  and  the  gradual 
increase  of  the  compressive  loading  area  in  the  core 
are  one  of  the  unique  features  in  the  deformation  and 
failure  behavior  of  the  deployed  foam-filled  composite. 
The  direct  consequence  of  the  extension  of  the 
damaged  and  compressive  zones  causes  the  simultan¬ 
eous  shift  of  the  neutral  axis  y*  and  the  damage 
boundary  yb.  Thus,  y*  shifts  from  its  initial  state  y* 

=  0  to  y*  <  0  as  yb  extends  from  yb  -  ~  to  yb  >  ~ 

when  o(f)  >  ob.  Both  y*  and  yb  move  downward 
during  the  extension  of  the  buckling  and  compression 
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zones.  The  amounts  of  their  extension  depend  on  E(s), 
and  G(f\  and  on  the  nonlinear  deformation  of 
the  core  material  which  can  be  determined  exactly 
from  the  present  analysis.  Figures  19,  21,  25,  30, 
36,  38,  44,  and  48  g've  quantitative  information  of 
the  distances  for  both  y*  and  yb  during  the  entire 
loading  history.  In  general,  it  has  been  obser  d  that 
the  extension  of  yb  is  much  larger  than  tha-  of  y*. 
For  example,  the  shift  of  yb  in  a  Kevlar  49  fabric/ 
4  lb/cu  ft  (64  k/m3)  polyurethane  foam  has  a  value 
about  half  of  the  diameter  of  the  beam,  while  y* 
shifts  only  less  than  one  tenth  of  the  cylinder  diameter, 
when  P  =  2000  lb  (909  kg). 

Additional  Remarks 

While  the  present  analyses  provide  accurate  means 
of  predicting  nonlinear  deformation  and  local  buckling 
of  composite  cylinders  subjected  to  four-point  bend¬ 
ing,  several  material  and  geometric  complexities  have 
not  been  considered.  These  may  include  the  anisotropy 
of  the  foams  and  fabrics,  the  heterogeneity  of  the  core 
material  and  the  fabric/urethane  composite  skin,  and 
the  geometric  nonlinearity  caused  by  large  strains 
developed  in  the  core.  Excessive  creep  of  the  foamed 
core  under  sustained  loading  has  also  been  noticed,  but 
is  not  included  in  the  present  study.  Effects  of  these 
complexities  on  deformation  and  failure  behavior  of 
the  composites  have  not  yet  been  fully  explored. 
Further  investigation  in  these  areas  should  be  of 
significant  assistance  in  a  better  understanding  of  the 
fundamental  behavior  of  the  material  and  will  enhance 
the  current  capacity  of  analyzing  and  predicting  more 
complicated  foam-filled  composite  structure  problems. 

6  CONCLUSIONS  At  ;D 
RECOMMENDAT'ONS 


Based  on  the  analytical  formulation  developed  in 
Chapter  3  and  the  results  obtained  in  Chapter  5,  the 
following  conclusions  may  be  drawn: 

1.  An  efficient  and  accurate  method  of  analysis, 
formulated  by  the  nonlinear  finite  element  method 
and  local  elastic  stability  theory,  has  been  developed 
successfully  for  studying  the  high-strength  fabric/ 
foam  composite. 

2.  Deformation  of  the  cylindrical  composite  struc¬ 
tural  element  and  progressive  local  failure  behavior 
of  the  load-bearing  skin  in  the  composites  with  various 
high-strength  fabric  materials  and  core  foams  can  be 


predicted  accurately  using  the  present  method  of 
analysis.  These  results  have  been  confirmed  by  experi¬ 
mental  data. 

3.  Transverse  shear  effects,  usually  considered  to  be 
negligibly  small  in  conventional  structural  members, 
are  shown  to  play  an  extremely  important  role  in  (he 
overall  deformation  of  the  composite  and  the  locitl 
failure  of  the  skin.  They  should  be  incorporated  intti 
the  formulation  to  provide  an  accurate  analysis  and 
prediction. 

4.  Local  buckling  (or  wrinkling)  of  the  load-bearing 
composite  skin  is  found  to  occur  at  a  very  low  loading 
level  and  to  initiate  nonlinear  deformation  of  the 
composite  in  the  present  four-point  bending  of  the 
composite  cylinder. 

5.  The  nonlinear  constitutive  relationship  of  the 
com  material  has  been  observed  to  significantly  affect 
the  deformation  and  local  failure  of  the  composite. 
This  is  especially  true  at  the  later  stage  of  composite 
deformation. 

6.  Progressive  extension  of  the  failure  zone  in  the 
skin  causes  continuous  redistribution  of  internal 
stresses  in  the  composite.  Direct  consequences  of  this 
progressive  failure  are  manifested  by  the  reduction  of 
the  load-bearing  cross-sectional  area  of  the  skin  and  the 
progressive  extension  of  the  compressive  domain  in 
the  core,  which,  in  many  cases,  are  so  significant  that 
the  neutral  plane  is  shifted  almost  to  the  lower  skin. 

Based  on  the  experimentally  observed  and  analyti¬ 
cally  predicted  performance,  it  is  recommended  that: 

1 .  The  high-strength  fabric/foam  composite  system 
be  used  as  a  load-bearing  element  in  a  certain  class  of 
structures. 

2.  The  currently  developed  analytical  method  for 
predicting  the  nonlinear  deformation  and  progressive 
failure  of  the  composite  be  used  in  the  analysis  and 
design  of  the  composite  structural  elements. 

3.  The  foam-fabric  composite  material  identified 
in  this  study  be  investigated  further  to  resolve  such 
features  as  ideal  section  design,  foam  density,  and 
fabric  design.  The  composite  appears  to  offer  an 
excellent  means  of  constructing  low  design  load 
structures  but  further  study  is  needed  in  joint  design, 
fastening  performance,  and  component  processing 
techniques. 
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Figure  18.  Deflection  of  central  plane  (at  z  =  0)  in  deployable  composite  cylinder 
subjected  to  continuously  increasing  load  (Kevlar  49  fabric/urethane  skin 
and  B-F  foam ;  d  =  12  in.). 


Figure  20.  Distribution  of  at  in  the  core  of  a 
composite  cross  section  at  z  =  0,  P 
=  1200  lb  (Kevlar  49  fabric/urethane 
skin  and  B-F  foam;d  «  12  in.). 


Figure  21.  Distribution  of  oz  in  the  core  of  a 
composite  cross  section  at  z  =  0, 
P=  1600  lb  (725  kg)  (Kevlar  49 
fabric/urethane  skin  and  B-F  foam; 


d  =  12  in.  [305  mm]). 


Figure  22.  Distribution  of  oz  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =  2000  lb  (907  kg)  (Kevlar  49 
fabric/urethane  skin  and  B-F  foam 
d  =  12  in.  [305  mm]). 


Figure  23.  Distribution  of  oz  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =2400  lb  (  ’8  kg)  (Kevlar  49 

fabric/urethane  skin  and  B-F  foam; 
d  =  12  in.  [305  mm]). 


Figure  26.  Distribution  of  oz  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =  500  lb  (227  kg)  (Kevlar  49 
fabric/urethane  skin  and  A-F  foam  ; 
d  =  12  in.  [305  mm]). 
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Figure  27.  Distribution  of  az  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =  800  lb  (362  kg)  (Kevlar  49 
fabric/urethane  skin  and  A-F  foam; 
d  =  12  in.  [305  mm]). 


Figure  28.  Distribution  of  oz  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =  1000  lb  (453  kg)  (Kevlar  49 
fabric/urethane  skin  and  A-F  foam ; 
d  -  12  in.  [305  mm]). 


sEses*  . 


Figure  31.  Distribution  of  az  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =  800  lb  (362  kg)  (Kevlar  29 
fabric/urcthane  skin  and  B-F  foam; 
d  =  12  in.  [305  mm]). 
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Figure  32.  Distribution  of  at  in  the  core  of  a 
composite  cross  section  at  z  -  0, 

P  =  1200  lb  (544  kg)  (Kevlar  29 
fabric/urethane  skin  and  B-F  foam; 
d  =  12  in.  [305  mm]). 


Figure  33.  Distribution  of  oz  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =  1600  lb  (725  kg)  (Kevlar  29 
fabric/urethane  skin  and  B-F  foam ; 
d  =  1 2  in.  [305  mm]). 


Figure  34.  Distribution  of  oz  in  the  core  of  a 
composite  cross  section  at  z  =  0, 

P  =  2000  lb  (907  kg)  (Kevlar  29 
•  fabric/urethane  skin  and  B-F  foam; 
d  =  1?  in.  [305  mm]). 
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Figure  35.  Deflection  of  central  plane  (at  z  =  0)  in  deployable  composite  cylinder 

subjected  to  continuously  increasing  load  (Kevlar  29  fabric/urethane  skin 
and  A-F  foam;  d  =  12  in.  [305  mm] ). 
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Figure  36.  Extension  of  failure  zone  in  skin  and  shift 
of  neutral  axis  in  a  composite  section  at 
l  =  0  (Kevlar  29  fabvic/urethane  skin  and 
A-F  foanv.d  =  12  in.  |305  mini). 
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Figure  37.  Deflection  of  central  plane  (at  z  =  0)  in  deployable  composite  cylinder 

subjected  to  continuously  increasing  load  (Kevlar  49  fabric/urethane  skin 
and  B-F  foam;  d  =  6  in.  [1S2  mm]). 


Figure  38.  Extension  of  failure  zone  in  skin  and  shift  of  neutral  axis  in  a  composite 
section  at  z  =  0  (Kevlar  49  fabric/urethane  skin  and  B-F  foam;  d  =  6  in. 
(152  mm]). 
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Figure  43.  Deflection  of  central  plane  (at  i  -  0)  in  deployable  composite  cylinder 

subjected  to  continuously  increasing  load  (Kevlar  29  fabric/urethanc  skin 
and  B-F  foam;  d  *  6  in.  [1 52  mm)). 
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Figure  44.  Extension  of  failure  zone  in  skin  and  shift  of  neutral  axis  in  a  composite 

cross  section  at  z  =  0  (Kevlar  29  fabric/urethane  skin  and  B-F  foam;  d  =  6  in. 
(152  mm)). 
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Figure  47.  election  of  central  plane  (at  z  =  0)  in  deployable  composite  cylinder 

subjected  to  continuously  increasing  load  (Kevlar  29  fabric/urethane  skin 
and  A-F  foam;  d  =  6  in.  1 152  mm]). 


Figure  48.  Extension  of  failure  zone  in  skin  and  shift  of  neutral  axis  in  a  composite 

cross  section  at  z  =  0  (Kevlar  29  fabric/urethane  skin  and  A-F  foam;  d  =  6  in. 
{152  mm]). 
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APPENDIX: 

ELEMENT  STIFFNESS  MATRIX 
k  OF  THE  8EAM 


i 

* 


k,  |  =C,  i/S2 

k12  =-(12|-6)C ,  2/C3 

k13--(6£-4)Cl2/8J 

k,4=-C  12/k2 

k,s  =-C,,/KJ 

k16=(12$-6)C12/«3 

k17  =-(6§-  2)C12/C2 

k18  =  C  j  2 /ft2 

k2J  =(12?-6)2  C22/84 

k23  =(12?-6)(6$-4)Cj2/K3 

k24  =  (12|  -  6)  C22/fi3 

k2S  =(12S  -  6)C12/e3 

k26  -  -  (6-  12|)2  C22/84 

k2  7  =  ( 1 2|  -  6)(6£  -  2)  C2  2  /83 

k28  =- (12^-6)  c22/e3 

k33=(6|-4)2  C2I/e2 
k34=(6$-4)C22/K2 
k3S  ::  (6|  -  4)  C 1 2 /C2 


k36 

=  (66  -  4)(6  -  1 2«)  C2  2  /C3 

k3? 

=  (6*-2)(6$-4)C22/R2 

k38 

=  -(66-4)C22/C2 

k44 

=  C22/82+(l-$)’C33 

k4  s 

=  c12/e2 

k46 

=  (6  -  12£)  C22/83 

k47 

=  (6$-2)C  22/82 

k48 

=  -c22/e2  +c33  |(i  -?) 

ks  s 

=  Cn/CJ 

kj6 

=  (6  -  12{)  C]  j/83 

ks7 

=  (6{-2)C12/82 

kss 

=  -C12/22 

k66 

=  (6-12{)2  c22/e4 

k6  7 

=  (6£  -  2)(6  -  12|)  C22/83 

k6s 

=  -  (6£  -  12)  C22/83 

k77 

=  (6£-2)2  C22/f5 

k3a 

=  -  (6£  -  2)  C22/82 

kss 

=  c22/e2  +  £2  c33 
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